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SUMMARY 



An analysis has oeea made of a Gonsidera"ble amount of 
data for turlDulent "boundary layers along wings and "bodies 
of various shapes in order to determine the fundamental 
varialJles that control the development of turbulent bound- 
ary layer It was found that the type of velocity dis- 
tribution in the boundary layer could be expressed in 
terms of a single parc^meterw This parameter was. chosen 
as the ratio of the displacement thickness to the momentum' 
thickness of the boundary layer. The variables that con- 
trol the development of the- turbulent boundary layer appar- 
ently are (l) the ratio of the nondimensl onal pressure 
gradient, expressed in terms of the local dynamic pressure 
outside the boundary layer and boundary-layer thickness, to 
the local skin- friction coefficient and (2) the shape of 
the boundary layer, Ax}. empirical equation has been devel- 
oped in terms of these variables that, when used with the 
momentum equation and the skin-friction relation, makes it 
possible to trace the development of the turbulent bound- 
ary layer to the separation point. 



INTRODUCTION 



A good measure of the understanding of the general 
problem of the flow of a real fluid about a body of arbi- 
trary shape can be taken as the degree of approximation 
with which the aerodynamic characteristics of the body can 
be calculated on the basis of existing knowledge. The 
flov/ in regions removed from the v/ake and from, the surface 
of a body obeys very closely the laws for perfect fric« 
tionless fluids at all reasonably large Reynolds numbers. 
This type of flow is well understood, although the detailed 
computations may .be difficult in some cases. Departure of 
the flow of real fluids from that of the frictionless fluid- 
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is caused almost entirely "by the failure of the idealized 
fluid to reproduce the actual flow conditions at the sur- 
face of the hody* 

Because the theory of perfect fluids gives zero drag 
for all hodies and gives no information concerning con- 
dibions that lead to separation of the flow from the sur- 
face^ these phenomena must be almost entirely associated 
\^rith the 'behavior of the flov; at the surface or, in other 
wordSj with the "boundary layer; that is, all drag, with 
the exception of induced drag and drag due to shock waves, 
all cases of flow separation ouch as occur at maximum lift 
or at high aileron deflections, and all Reynolds number ef- 
fect b are entirely dependent .on the "behavior of the "bound- 
ary layer. 

Because the flow of ideal fluids is well understood, 
the proDiem of calculating the actual flov/ resolves it- 
self into a study of boundary layers. . There are three 
general types of boundary layer; laiainar, transitional, 
and turbulent. In most cases^ the important laminar 
boxmdary'-layer characteristics - thickness, skin friction, 
and point of separation - can be estimated with sufficient 
accuracy from methods described in references 1 and 2. If, 
in unusual cases? mors detailed information Is required, 
one of the methods of calculation given in reference 3 may 
be used« 

Experimsut s in 1 o w-^tur bulence v/ind tunnels and in 
flight have shown that extensive laminar ].ayers may be 
maintained in the presence of a favorable pressure gradi- 
ent. Although little is known about methods of determin- 
ing the position of the transition point, transition must 
occur either associated with laminar separation or at some 
position upstream of the separation point. The position 
of the transition point varies widely within this region 
because of changes in stream turbulence, surface condition, 
and Eeynolds number o 

A considerable amount of progress has been made in 
determining turbulent skin friction in pipes and along 
flat plates with zero pressure gradient. As a result of 
the development of the "mixing-length" theory, notably by 
Prandtl and von Karm^n, • reliable rules that should be ap- 
plicable throughout an extremely wide range of B>eynolds 
numbers have been found for calculating the skin friction 
and velocity distribution in the boundary layer along 
flat plates and in pipes. 
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In order to find the effect of varying pressures 
on the "boundary-layer characteristics, von Kirmin applied 
the momentum theorem to the "boundary layer and derived 
the so-called '^momentum equation", that gives the rate of 
thickening of the "boundary layer if the type of velocity 
distrihution within the "boundary layer, the external pres- 
sure distrihutioh, and the skin friction are known. It 
was found in a number of cases that, if. the shape of .the 
boundary layer and the skin friction were assumed to "be 
the same as in pipes, good agreement was obtained betv;een 
the calculated and experimental boundary-layer thickness 
and skin friction. This procedure, however, fails to 
•give any information concerning the changes in. boundary- 
layer, characteristics that lead to separation of the flow 
from the surface. In addition to the momentum equat i on , 
a relation is needed betv/een the shape of the velocity 
distribution in the boundary layer, the skin friction, and 
the pressure distribution. 

Several attempts have beisn '.made to find such a rela- 
tion., ffrom experiments in converging and diverging chan- 
nels, Flkuradse, using water (reference 4), and Donch, 
using air (reference 5), found that, in cases in which the 
boundary layer along, the walls met in the center of the 
channel, the velocity distribution across the channel -was • 

a function of a vR, where a- is the angle of diver- 
gence of the channel and R is the Reynolds number based 
on the channel width and average velocity. 3uri, whose 
work was discussed by Prandtl in reference 6, used the 
results of Kikuradse and Donch in an attempt to calculate 
the general behavior of turbulent boundary layers. Buri 
assumed that the shape of the boundary layer was always 
given by the -parameter of Nikuradse and D'onch. One of the 
weaknesses in Buri»s calculations was the assumption t^hat 
the shape of the velocity distribution depended only on 
the value of this parameter and was independent of the 
previous history of the boundary layer, "it may- be pointed 
out that the parameter of Hikuradse and Donch is not essen- 
tially a shape par^ameter but represents a function of the 
pressure gradient and skin friction which, is assumed to 
determine the shape of the velocity distribution. 

&ruschwitz (reference 7), using a parameter dependent 
upon the type of velocity di stribution. in the boundary 
layer, found a relation between the pressure gradient and 
this parameter. This relation, with the momentum equation, 
was sufficient to determine the d-evelopment of the turbu- 
lent boundary layer along a surface. Although Gruschwitz 
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obtained good agreement with experiment for .the data pre- 
sented in his paper, other, investigators v/ho have tried 
to use the method reported .poor agreement for cases in 
-which the turbulent "boundary layer separated from the sur- 
face. Peters (reference 8). conducted an investigation 
for the specific purpose of. testing the ' Gruscliwit z method 
of calculation and , concluded that the Gruschwitz method 
cannot "be tised to determine the location of the separation 
point nor even^ in many cases, tO; predict whether sepa- 
ration will occur at all.. 

The purpose of the present investigation is. to deter- 
mine the important variables that control the behavior^ of 
turbulent boundary layers and to develop general relations 
in terms of 'these variables that describe the boundary- 
layer motion. The fundamental variables must, of. course, 
be expressed nondimensionally in terms of local boundary- 
layer quantities. 

The quantities at a given station along- a surface 

that v/ere felt to have the niost important effect on the 

further development of the boundary layer are the follow- 
ingt 

(1) Shape of the boundary-layer ' prof ile 

(2) Rate of change along the surface of the dynamic 

pressure outside the boundary layer 

(3) Skin friction 

Data from various published sources and from tests in the 
HACA .two-dimensional low-turbulence tunnel made specifi- 
cally for the present investigation v;ere analyzed in terms 
of the foregoing variables in order to find the needed 
general relation for the rate of change of shape of the 
boundary- layer profiles * 

■SYMBOLS 



a angle of divergence of channel; also, angle of attack 

|l Reynolds number 

D minor axis of ellipse 

u velocity within boundary layer 
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\Jq free-stream velocity 

X distance along surface 

c chord 

y distance perpendicular to surface 

e momentum thickness | / ^ (}' ^} ^^1 

J 0 

U velocity outside "boundary layer 

skin friction per unit length 
q dynamic pressure outside "boundary layer 
E .shape parameter 



6 displacement thickness 

p density 

]}» viscosity 

6 "boundary-layer thickness 
V ' kinematic viscosity 

C constant ' 

s equivalent length of flat plate "before pressure re- 
covery "begins 

Id width of channel 

hj. total pressure at distance 6 from surface 

h^ free-stream total pressure 

T] shape parameter 

Uq velocity at distance 6 from surface 



r /■"Iv 1 
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initial value 




<1 
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constant 
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initial value 


of 


■ e 




initial value 


of. 
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2.557 logg 4. 
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2.557 . 1 
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SXPEHIIOHTAL ' DATA « 



The data used in the analS'sis were collected from the 
available literature and from tests performed in the NAOA 
two-dimensional 1 ow- turtulence tunnel. The following ta- 
hie shov/s the data used in the present analysis; 



Model 


fieynolds 
nunher, E 


Angle of 
attack 
(deg) 


Reference 


NAG A 66,2-216 


0.9 X 10^ 

1.5 

2.2 

2.5 

2.6 


10.1 
10.1 
10 .1 
10.1 

10,1 


9 


STACA 65(216)-222 
(appr OS. ) 


.92 
1.51 
2.67 

.92 

1.51 
2.67 


8.1 

• 8.1 
8.1 
10,1 
10.1 
10.1 


Present 
report 


HACA nose- 
opening 

airfoil shape 13 


1.46 
• 2.39 
4.18 


9 .1 
9.1 
9.1 


Present 
report 


Channel 




Eun 2 


7 


JIACA 0012 


7,6 


0- 


10 


Symmetrical 
airfoil 


3 .8 (approx . ). 


9 


3 


Slliptic cylinder 


.118 


0 


11 



The paper "by &ruschwit z ,. ( ref ef ence 7) contained" data ph a 
v;lng and on channel wails /©"btaihed from te st s at G-ottingen. 
The' data from run 2 on.the 6hannel wall in the forjfn of 
boundary-layer velocity prof ile s , pressure di st rihut i on , 
and curves of the momentum. thickness and the shape-param- 
^■ter plotted against distance along the plate were used. 
Prom Peters' paper (reference 8), data were taken 'in- the 
form of curves of the mpmentum thickness, shape [parameter, 
and pressure coefficient against the position along 'th.e;^'* 
air foil chord p . These tests were made in order to-cn'eck- 
the method of calculation proposed by. Gruschwi t 2 ♦ The 
Reynolds numher of the test vras not given explicitly and 
was judged, to he slightly less than 4^000, 000. The re- 
sults from the tests at an angle of attack of 9 were used 
"because the data for this angle of attack were presented 
in a convenient form and "because separation of. the flow 
had taken place at the rear of the wing, 

A few points v;ere obtained from the houndar3r-layer 
velocity profiles and pressure distribution contained in 
reference 11, Th*e* data v;ere obtained from a test of "an 
elliptic cylinder at an angle of attack of 0^ and a *• 
Heynolds number of 118^000 based on the minor axis D of 
the ellipse., The ratio of major to minor axis v;as 3.95, 
making the Seynolds number based on the major axis eqiial . 
to 349,000. 

The data on turbulent boundary layers involving sep- 
aration obtained from the KACA tv/o-diiaensional lov/- 
turbulence tunnel T^ero collected from previous tests and • 
from tests performed specifically for t;he present investi- 
gation. 

A fevr turbulent boundary-layer profiles that were not 
clc.se . to separation were obtained from previcusiy published 
data on the iTACA 0012 airfoil (reference lO) . These points 
were used in the anc^ly^^is mainly because the boundary layer 
v/as far f r cm " separat i cii and the use of these point g helped 
to ^i'^^e a better di Btr .1 Vat 1 on of data^ The data obtained 
on the >TAGA 66,2-216 airfoil are given in .reference' 9, Or- 
dinr.tas for this section can be found by methods described 
in reference 12, 

Tests of a thick airfoil, the KAGA 65(216)-222 ( ap- 
prox.. ) v/ere made at three Heynolds numbers at tv^o angles 
of attack. The methods of obtaining tjie ,data were the 
same as those described . in ref e>enQe;:9 . / All ' tfie^- test^^ of 
this airfoil involved turbulent separation^' The data 
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from these tests are presented in the form of "boundary- 
layer velocity profiles for a number of stations along 
the chord (figs. 1 to 4). The pressure distributions are 
given^ in figure 5» The region of turbulent separation is 
indicated in the pressure* di stribut ion as the flat region 
at the rear of the airfoil. The chord of the airfoil was 
24 inches and the airfoil v/as finished as described in, 
reference 9o The ordinates can be derived by methods ex- " 
plained in reference 12. The finish \7as free from all 
surface imperfections that could be felt by hand but had 
a strip of carborundum-covered cellulose "Scotch" tape.l 
inch wide on the upper • surf ace near the leading edge. 

The other model tested in the UACA two-dimensional 
low-turbulonce tunnel for data to be used in the present 
analysis was the NACA nose-opening airfoil shape 13. The 
model had a chord of 36 inches and was finished in the 
same manner as the NACA 66,2-215 (reference 9). The sec- 
tion ordinates for the ITACA nose-opening airfoil shape 13 
are given in reference 13, The wing was tested at an angle 
of attack of Qd^ at three Reynolds numbers. The turbu- 
lent separation obtained in this test was not so marked as 
that obtained in the tests of the HACA 56,2-216 and xTACA 
65(216)-222 (approx.) airfoils, although tufts placed at 
the rear of the v/ing on the upper surface indi cat ed ' sepa- 
ration. The data from these tests are presented in fig- 
ures 6 and 7 in the form of boundary--layer velocity pro- 
files for a number of stations along the chord. The'pres- 
sure distributions are given in figure 8. The beginning 
of separation is indicated in the pres55Ure distributions 
by the flattening of the curves at the rear of the airfoil. 
The small flat region in the pressure distribution at the 
nosG of the airfoil is an indication of laminar separa- 
tion. The boundary-layer velocity profiles for the region 
at the nose of the airfoil are shown in figure 6. The pe- 
culiar shape of the velocity distribution for some of the 
stations,, particularly in the curves that show increasing 
velocity with approach tov;ard the wall, is probably caused 
by spanv;is6 flov;s over the airfoil. The boundary-layer 
thicknesses obtained in these tests were muph larger than 
those usually obtained for airfoils of 35-inch chord. 



ANALYSIS 



The equation that gives the rate of change of the mo- 
mentum defect in a boundary layer, originally derived by 
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von K^rman, may "be written in the follo^^^ing form for two- 
dimensional flow: 

dx \ 2 / q dx 2q 

where 

CO 

8 momentum thickness \ j ^ " 

u velocity i^ithin lioundary layer 
U velocity outside boundary layer 
y distance perpendicular to surface 

skin friction per unit length 
q dynamic pressure outside "boundary layer 
X distance along surface 
H shape parameter 

The shape parameter E is defined as the ratio 6*/6 

/"i CO 

where 6* j (l ^ dy. The difference between the ac- 

0 

tual flow of momentuiu in the boundary layer and that of * 
the same quantity of fluid flowing with velocity U is " 
pU e. From this relation, the length 6 is given the 
name momentum thickness..-. The length 6*, called the dis- 
placement thickness, is the amount by which the stroam- 
lines just outside the boundary layer are displaced because 
of the reduction of ^velocity within the boundary leyer. 
Because 6 depends "on the second power of the velociity 
distribution, whereas 5* depends on only the first power, 
the ratio S*/8 depends on the manner in which u/U var- 
ies with y - that is, upon the shape of the boundary-layer 
profile » 

The momentum equation in the form just- ^iven contains 
only local boundary-layer quantities. The local skin- 
friction coefficient Is To/q. The nondimensional i^ressure 
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gradient is - -^•S' , where 9 is the unit of length and 

O CX ISC 

q is the unit of dynamic pressure. 

Although it has heen shown that the shape of the 
houndary-layer profile determines H, the converse cannot 
"be proved from mathematical considerations alone. If H 
does actually determine the sliape of the boundary-layer 
profile, then all points of u/'u plotted against H at a 
constant value of y/8 should fall on a single curve, A 
collection of such curves for various values of y/9 is 
shov/n in figure 9o The data presented in figure 9 repre- 
sent the collection of all the boundary-layer .profiles 
that enter into the analysis* figure 9 shows, that u/U is 
a function of H alone for a given value of y/Q. 5?his • 
conclusion is important because ■ it means that turbiilent— 
boundary-layer profiles fori?, a single— par ame ter family of 
curves o The complete velocity d is tr ibu.t i on in the bound- 
ary layer is known when 9 and H have been determined. 
Figure 10, which is a. cross plot of figure 9, gives turbu- 
lent-boundary-layer velocity profiles corresponding to 
various values of- As the separation point is approached, 

the value of H increases. Because the turbulent separa- 
tion point usually is not very well defined,, it is not pos- ' 
siblo to give an exact value ox H corresponding to sepa- 
ration. The value of H usually varies so rapidly near 
the separation point, however, that it is not necessary to 
fi>: accurately the value of H corresponding to separa- 
tion. Separation has not been observed for a valxie of H 
less than 1.8 and appears definitely to have occurred for 
a value, of H of 2,6e Gruschwits's criterion for immi- 
nent separation ia eqtuivalent to a value of H of 1.85. 

The fact that the typo of ve 1 oc i ty " di s tr ibut i on in 
the boundary layer can be -given in terms of a single pa- 
rameter greatly simplifies the study of turbulent boundary 
layers* It is now necessary to determine only the manner 
in v/hich this' parameter varies a.long the surface as a 
function of the external forces a.cting on the boimdary 
layer » 

The external f or co s - ac t ing on the bound?iry layer at 
any point are the pressu.re gradient, expressed nondimen- 

sionally as - ai^d the skin friction, expressed non- 

a dx 

dimensionally as The ass-ximption is made that the 

rate of change of H rather than PI itself is related to 
the local forces. This assumption is desirable in order 
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that the hounda.ry-layer conditions downstream from a point 
shall he definitely connected with the conditions upstream 
of the point; that is, a sudden change in the pressure 
gradient should not produce a discontinuity in the type of 
velocity distrihution in the "boundary layer. Prandtl (ref- 
erence 6) has pointed out this difficulty both in the Pohl- 
hausen theory for laminar, "boundary layers and in Buri*s 
method for calculating turbulent "boundary layers. Expressed^ 
in nondimensional form, the rate of change of H is given 

as e 1^- 

dx 

In the early stage of the analysis, the experimental 

data were plotted in the form of ' G 4^ against ^ 

dx q dx 

Fair correlation was o"btained for a limited amount of data. 
As the analysis \tf-as extended to include more data, syste- 
matic variations with Reynolds number were noticed. When 

9 — was plotted against — -f^ ~f the consistent varia- 
ax q clx 

tion with Reynolds number was eliminated^ The skin fric- 
tion ^^ras tentatively assumed to be given by the Squire 
and Young formula (reference 14) 

1^ = 5,890 log^^ (4.075 RQ)j 

where 

This formula was chosen because of the good agreement ob- 
tained between the experimental drag coefficients^ for air- 
foils and those calculated by the Squire and Young method. 
It was felt that the local skin-friction coefficient s thus 
determined were probably more accurate in m.ost cases than 
those determined directly from the boundary-layer surveys. 

4. . ^ 6 dq 2a 

ihe quantity ^ ^ quant ity . differing only 

by a constant factor, has frequently been used as a param- 
eter for boundary-layer phenomena. For example, the Pohl~ 

. ft 2 J TT 

hausen parameter for laminar boundary layers — — » where 

V dx 

6 is the boundary-layer thickness and v the kinematic 
viscosity, can be shoxvn to be equivalent to - — For 

a laminar boundary layer, t. = ^J^ (^] and is pro- 
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portional to T. It thus follows that — is -oroiDort i on- 

T T Ui 

al to -Tr- Substituting for in the Pohlhausen 

dU U dq 

parameter and replacing ~ hy its equivalent ~ 

gives a! quantity proportional to the Pohlhausen parameter 
g 6 dq 

q dx ~ 
6 dq 2q 



C 3- -r^, which in turn is equal to some constant tiaes 

q dx 



By a Similar process of reasoning, the parameter 
q dx 

— — , vhere s is the equivalent length of the flat plate 
U dx . . ' 

"before pressure recovery is "begun, which determines the 
amount of pressure that can he recovered in a laminar layer 
with a straight-line velocity gradient as given in refer- 

6 d Q 2 0 

ence 1, can "be shown to "be proportional to — — - • 

q To 

Nikuradse (reference 4) found that his results for turbu- 
lent flov;, which give the velocity distribution across 
diverging or converging channels , agreed xvith similar meas- 
urements made by Donch (reference 5) v;hen a 

same value. For a given type of velocity distribution 

b do 

across the channel, ^ at the center is proportional 

q dx 

to the angle of divergence, where b is the v;idth of the 
channel* Within the range of Reynolds numbers covered by 
the investigations of Niicuradso and Donch, the skin-fric- 

tion coefficient — ^ at the wall was inversely proportion- 
al to the V H. Tor a given type of velocity distribution, 

aVs is therefore praporiional to — 4^ — ~ • It niay be 
^ ' ; ■ q dx Tq 

pointed out that, although in 3uri»s theory of turbulent 
boundary layers and Pohlhausen' s theory- of laminar boundary 

6 dq 2o 

layers — — - x^ras assumed to determine the type of v€j- 
q dx 

-locity distribution in the boundary layer, in the present 

analysis it is assumed that 1 affects only the 

q dx T>Q 

rate of change of /the .type of velocity distribution. 

i' _ • • 

It seemed highly probable that the rate of change of 
H shoul.d depend.'not only on the ratio of the pressure. 
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gradient to the skin friction "but also on tho value of H 

itsolfii Plots wero thorefore made of e -r- against 1 ?S £i 

dx ^ q ax Tq 

at Ool intervals of H for all the data entering into the 

present analysis. These plots are given in figure ilo 

Although the points show considerable scatter^ definite 

trends for the variation of 9 ~ with both 1 £i and 

q dx Tq 

H are observableo It may be pointed out ^ however 5 that 
dH f\ 

^oth ^ and ^ were the slopes of experimentally deter- 
mined curvesg Largo scatter of the data therefore is to bo 
expected^ Tho large scatter of the points in figure 11, 
consequent ly 3 does not necessarily indicate any serious in— 
adoquacios in the present analysis « Prom a study of the 

available data, it was found that the variation of 6 ^ 
9 do 2x3 

with - 4 and H could be fairly well reprosonted by 

q dX Tq 

the equation 



^ dH 4 e 680 -a . 97 5 ) 

8 *T — = e 
dx 



^ 0 M |a . 2.035 (H - 1.286) 
_ q dx 



The exponential foriia of the factor multiplying the 
second member of this equation v/as chosen because the data 
for high values of H, although not very complete, never- 
theless indicated that ~ was large. It may be noted 
that 1.286 is the value . of H for the ■i'-power distribution 

of velocity in the boundary layer. It is seen from the 
foregoing equation that no change in H is indicated for 

the case in which ^ = 0 and H has the value 1.286o 

dH 

The degree to which the equation for 6 approximates 

the experimental data from which it v/as derived may be 
seen in figure 11. Each of the straight lines in figure 

11 vas obtained from the equation for 6 ^ by giving H 
successive values of 1.35, 1.45,- 1.55, and 1.65. The slope 

of these lines is given by the factor - eeo{E ^2 . 97 5) ^ 
and the intercept for Q = 0 is given by the term 
«2.035 (H - 1.286) • ^ 
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COMPARISOH OP PRESENT ANALYSIS WITH 
GEUSCHWITZ ANALYSIS 



In reference 7, Gruschwitz analyzes the "bshavior of 
turbulent "boundary layers in the following manner: 



e dh^ 
dx 



= I Cn, He) 



where 
and 



shape parameter 



in the saiie sense as li 



velocity at y ^'Q 



total pressure at distance 6 from surface {g^ in 
reference 7) 



In writing this relation, it is implicitly as^med that 



e 



— ~ IS independent of 



Because of this assump- 



dx - - - dx' 

tion, the G-ruschv/itz analysis is- suhject to the same crit« 
icism as are the Buri and Pohlhausen methods for making 
"boundary-layer calculations. Kow 



- h^ 



where 



pressure; hence, 



(go in reference 7) is the free-stream total 



e dh 



and therefore 



q dx 



dx 



„ drr 6 dq_ ^ 
dx 



q dx 



Except for the factor 



e dq 
q dx 



r, - F (Tl. Eq) 



- this relation, is similar in 



form to the relation found from the present investigation 
namely, 



dx . \q dx T,, J 
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The Gruschwitz relatioiii however, is very restricted in the 

form of the dependence of • 0 4^ on — 4^ in comparison 

dx q dx 

vrith the type of relation used in the present investigation. 

Furthermore, no variation of 6 with Ro is indicated 

dx 

in tho final equations given by Gruschwitz, 

In the Gruschwitz analysis, the arbitrary function 
v/hich was to he determined experimentally contained only 
one variable T) ; whereas, in the present analysis, the 

arbitrary function contains two variables, ^ ^ ^ and H. 

q dx 

Apart from the neglect of Eg, one reason for the failure 
of the Gruschwitz analysis is that a correlation of all 
turbulent-boundary-layer data in terms of a function con- 
taining only one independerit variable was not possible. 

METHOD OF CALCUIiATIOH 



For calculating the characteristics of turbulent bound- 
ary layers, the following information is required: The in- 
itial values of 9 and H, the pressure distribution over 
the body, and the Reynolds number. The equations that ar<3 
used in making a computation are 



+ £t2 e dq ^ 
dx 2 q. dx " 2q 



2l ^ 



^0 



5.890 loe.^^ (4.075 Sg) 



T 

In order to reduce tho work of computation, -r^ was plot- 



ted against Eg in figure 12 and the factor '&bo(^-s. 975) 

^ jn 

13. 



in the equation for || was plotted against H in figure 



The momentuin equation and the equation for are 

simultaneous firet-order differential eauations that can be 
solved by a step-hy-step calculation. It is usually neces- 
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sary to use such a method althou.gh, for some particular 
casesj the equations may "be integrated directly. The meth- 
od of calculation is as follows: The values of the var- 
iables entering into the computation at the initial sta- 
tion are substituted in the momentum equation and the equa- 
tion for Values for M and ~ are thus obtained 
dx dx dx 

at the initial station. An increment of the length along 
the surface of the body x is then chosen and multiplied 

by — and — to give AG and AH, respect ively« 
dx dx 

These increments of B and H are added to the initial 
values and result in. values of 6 and H for the new 
value of X. The process is repeated until the desired 
result has been attained. Separation may be considered to 
have occurred when H rises to about 2*6. 

The choice of the increment of x is a matter for 
the judgment of the individual investigator. As a general 
rule, the increments of- x' should bo made small v/hen 

or — changes rapidly from one value of x to the 
dx dx 

next. In order to decrease the length of the calculation, 
the increments of x must be chosen as large as is com- 
patible with the accuracy desired. For the computations 
that v/ere made in order to check the method of calculation, 
the increment of x for one step was so chosen that 

dH 



dlE 



^ (dx j ^ 0.0025, where A is the change 

between two successive values of x and' Ax is the incre- 
ment of X. This criterion furnishes a measure of the 
maximum error that can be expected* in AH for one step of 
the computation. When the flow approaches separation, H 
usually increases very rapidly and, in such cases, the 
foregoing criterion may he disregarded without appreciable 
error in the position of the separation point. By disre-' 
garding the criterion v/hen the flov; is close to separa- 
tion, the length of the computation may be reduced. 

A sample calculation for the HACA 65(216)-222 (approx.) 
airfoil section at a = 10.1^ and at H = 2.67 x 10® is 
given in table I. 

If the question of separation is not involved and if 
the variation of H along the surface is not of interest, 
.reasonably accurate values of 6 may be -obtained by assum- 
ing a constant value of H and merely using the momentum 



I 
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equation together with the skin-friction relation to de- 
termine 6. This procedure is substantially the sase as 
that of reference 14 where a constant value of H of 1.4 
was chosen for calculating the profile drag of airfoil 
sections* 

TESTS 0? METHOD OF CALCULATION 



In order to obtain a general check of the method of 
calculation and to determine whether the scatter of the 
points in figure 11 v/as primarily due to the difficulty in 
obtaining the slopes of experimental curves or to serious 
inadequacies in the analysis, computations A^j-ero carried 
through, for eight cases. For all the se computat ions , the 
initial values of K and 6 were obtained from experi- 
mental data. 

A comparison between the calculated and experimental 
variations of H and 6 along the surface is shown in 
• figures 14 to 24, A comptit at i on made by the Orvischwitz 
method (from reference 8) is included in figure 17. For 
some of the cases, comparisons of the calculated and ex- 
perimental boundary-layer profiles at one or two positions 
are also presented. In general, the calculations are in 
good agreement with the experimental curves^ Ho systemat- 
ic differences were found between the calculated and ex- 
periment al . curve s of H. Although the agreement between 
the calculated and experimental curves of 6 (figs. 14 to 
24) is good in most cases, some consistent differences are 
apparent as the separation point is approached. In this 
region, many of the calculated values of 8 are less than 
the experimental values. One explanation of the discrep- 
ancy, of course, is that the Squire and Young skin-friction 
relation is in error in not indicating an increase in skin 
friction as the boundary-layer velocity profile approaches 
the shape for separation. This tendency, howev^er, is con- 
trary to the general impression that the skin friction 
should decrease as the separation point is approached. As 
the flow approaches separation, the fluctuations in the di- 
rection of flow increase. Such fluctuations make a pitot 
tube read velocities higher than the actual \^elocities. 
These fluctuation s are a large proportion of the mean flow 
close to the surface where reversed flo^-^ first begins. 
This behavior of a pitot tube may explain why the turbu- 
lent velocity profiles, v;hich are close to the separated 
st'ate, all have the characteristic hump at small values of 
y/c. The velocity profiles for lar.c;e values of H are 
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therefore in error in the region close to the surface. The 
error in profile shape affects 6 as well as H. Reading 

velocities too high for all points where ^< 0.5 makes 
the integral for 6, which is ■ ^1 ^ ^ C^}' large] 

than it actually is. It is to be expected, therefore, that 

rvalues of 6 determined "by pitot-tu"be readings should he, 
higher than the true values under conditions of unsteady 
flow. On the other hand, as separation is approached, the 
relatively greater velocity fluctuations near the surface 
may cause the skin friction to he higher than when condi- 
tions are far from separation; and the effect of the lower 
average velocities near the surface, such' as occur for 
higher values of • K, may thus oe compensated and possibly 
overbalanced. An attempt was made to- correlate the local 
skin-friction coefficient with PI, but no consistent re- 
sults v/ere obtained. Although there still' is considerable 
doubt concerning the true value of the skin-fri ct ion . coef- 
ficient for conditions approaching separation, it is in- 
teresting to note that the Squire and Young skin-friction 
relation was used through a range of Eg from 500 to 
48,000 and apparently gave good results for most of .the 
region covered by the turbulem; boundary lawyer. 

In making a calculations the initial condition of 
the boundary layer must be knovm* Computations which have 
been made do not indicate that the calculation for the 
curve is especially sensitive 
When the calculation is to be 

strong 



boundary layer is in a 

dx 



H- 

^0 the initial value of 9. 
made for a case in v/hich the 
adverse pressure gradient - 



is 



that 

than 2.035 (H - 1.28S) 
accurately determined. 



is of the same order or greater 



- the initial value of H must be 
This fact is obvious because the 
value of E determines how close the flov/ is to separa- 
tion. If the calculation is begun in a region where 



dx 



not very 

if i3. = 

dx 

gardless 



is very small or positive, the boundary layer is 

sensitive to the initial value of H. For example, 
0, H will eventtially have the value of 1,286 re- 
ef its initial value • 
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The "boundary-layer thickness is not particularly sen- 
sitive to the initial value of H. This effect is easily - 
verified from the form of the momentum equation. The ef- 
fect of changes in the initial value of 6 on the hound- 
ary-layer thickness further dovrnstream depends on the rel- 
ative magnitude of *5ie pressure . gradi- 
ent is large in comparison v;ith the skin friction, a 
change in the initial value of 8 v/ill produce a proper^ 
tional change in the subsequent values of 6 vrhereas, if 
the pressure gradient is small in comparison with the skin 
friction, a change in the initial value o.f 9 v;ill pro- 
duce a constant increment equal to the initial change • 

The initial value of 6 ordinarily may be taken the 
same as the value for the laminar boundary layer at the 
transition point. I:Iot enough is fcnovm about the mechan- 
ism of transition to be able to state in general what the 
initial values of H should be. If transition occurs 
in a zero or "favorable" pressure gradient, or if the 
joundary layer is sufficiently thin at the transition 
point, in accordance with the foregoing discussion, the 
position of the turbulent separation point will not be 
greatly affected by the choice of the initial value of ^H, 

PIIBSSUHE HECOVEEY WITHOUT CHANGE IN. 
BOUNDARY -LA YBE SHAPE 

dH 

The equation for 6 ^— indicates that, for each value 

bi il, such a value may be assigned to - ~ that 6 4^ 

q cLx ax 
= 0. Pressure may be recovered, therefore, v/ithout a 
change^ in boundary-layer shape if the correct pressure dis- 
tribution is used. The necessary pressure distribution 

can be obtained by using the equations for ^ and 

dx dx 

These equations can be integrated directly If H is as- 
sumed to have a constant value. When 4^ = 0. 

dx ' 

I I; = - ^'^^^ - ^-285) = Ci (1) 
where Cj is a constant. The momentum actuation then re- 
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duces to the form 



where 



dx 2q 



C = 1 + 2.035 (^ ^ 3 ^ ) (H - 1.286) 

T 

Sliminat ioil of -r^ hetv/eon equations (l) and (2) gives 

2q 

the relation 

• • 1 de ^ 6 dq 
C dx Ci q dx 



Upon integration, 



C 



«-o 



(3) 



vrhere Bq is the initial value of 9 and q^ the initial 
value of a» The relations hetween G and x and, con- 
sequently, "between q a,nd x are ohtained from integra- 
tion of equation (2). 
T 

For — , the Squire and Young skin- xr i ct ion relation 

2q 

v/ill he used 



2^1 [2.557 lege (4.075 Rq)V 

.^e-e.(t)(tr. 

where Eg is the initial value of Hp or 



^6 = ^6, 
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19. 

2q. 



, , , _, 

|2.557 lege 4.075 + • + l) log^ -^Jj 
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|3,557 l^loge 4.075 Rq^ + + logg 



9 



2 



1 

r 



The variables are separaljle. Let • 



= 2.557 (-^ + 1^ B = 2.557 logg 4.075 Rg^ 



ft X 

and, v;hen = 1, ~ 0 

^0 ^0 



The final equation then hecomes 



vjhere 



C = 1 + 2.035 ( — ^) (h - l,286j 

the use of equations (3) and (4), curves of ~ 
and S- against which indicate the inannei- in which 

^0 ^0 

presstire can he recovered without a change in "boundary-; 

layer shape, may he plotted, Figure 25 is a plot of £i 

C 

against ' 3, which shov;s the amount of pressure that can 
he recovered for a given change in 6 as a function of 
the shape profile to he maintained. The plot indicates 
that pressure can he recovered at the most rapid rate for 
a value of H fairly close to the value for separation; *' 
that is, ahout 2*3 . Because the flow with such a high 
value of K is apt to he unsteady, a good compromise he- 
tv/een steady flow and minimum increase in 6 with de- 
crease in q would seem to he a value of H of 1.7 or 1*8 
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OONCLUSIOITS 



An analysis of a considera/ble amount of data for tur- 
bulent "boundary layers collected from the availal>le lit-- 
erature and from tests in the KACA two-dimensional . lov/- 
turbulence tunnel Indicates the following conclusions: 

1. The shape of all . turbulent boundary-layer pro- 
files can be expressed as a function of a single param- 
eter • 

2. The variables that control the development of 

the turbulent -boundary layer apparently are (l) the ratio 
of the nondimensional pressure gradient, expressed in 
terms of the local dynamic pressure outside the boundary 
layer and boundary-layer thickness, to the local skin- 
friction coefficient and (2) the- shape of the boundary 
layer . 

3. An empirical equation has been developed in terms 
of these variables that, v/h&n used with the momentum equa- 
tion and the skin- f r 1 ct i on relation, makes it possible to 
trace the development of the turbulent boundary layer to 
the separation point. 

4. ITo systematic variation of the skin-friction co- 
efficient with the shape parameter was indicated by the 
data. 

5. Separation occurs for values of the shape param- 
eter greater than 1.8 and less than 2.6. 



Langley Memorial Aeronautical Laboratory. 

National Advisory Committee for Aeronautics, 
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Figure 16. -Comparison of experimental and calculated velocity profiles. 
(Data from reference 7, run -2.) 
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(a) Shape parameter H. 

Figure 17. - Variation of experimental and calcvdated values of ahape parameter H and 

momentum thickness e/c with x/c for symmetrical airfoil. R, 3.8 x 106; a go 
(Data from reference 8.) * > 7 . 
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Pigure 25* - Variation of experimental and calculated values of shape parameter H 
and momentxua thickness 6/c with x/c. Airfoil section^ NACA 63(21b)-222 (approx«)» 
R, 2,67 X 106; a, 8.10* 
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